Genuine collective flow from Lee- Yang zeroes 
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We propose to use Lee- Yang theory of phase transitions as a practical tool to analyze experimen- 
tally anisotropic flow in nucleus-nucleus collisions. We argue that this method is more reliable than 
any other method, and that it is the natural way to analyze collective effects. 
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Fifty years ago, Yang and Lee Q] showed that phase 
transitions can be characterized by the locations of the 
zeroes of the grand partition function in the complex 
plane. Since then, their theory has been extensively used, 
in particular, to study phase transitions in finite-size sys- 
tems, via numerical simulations |2j: in lattice calcula- 
tions, it has been applied to the electroweak || and QCD 
phase transitions 

In this Letter, we propose to apply Lee- Yang theory 
for the first time to the analysis of experimental data. 1 
More specifically, we show that it is the most natural 
way to study anisotropic flow in nucleus-nucleus colli- 
sions. Anisotropic flow is defined as a correlation between 
the azimuthal angle </) of an outgoing particle and the az- 
imuthal angle $ r of the impact parameter (see Fig. ^ 
$>r is also called the orientation of the reaction plane), 
which is best characterized by the Fourier coefficients of 
the single-particle distribution [fj: 



v n = (cosn(<j> - & R )) 



(1) 



In this expression, n is a positive integer and angular 
brackets denote an average over many particles belong- 
ing to some phase-space region, and over many collisions 
having approximately the same impact parameter. In 
particular, the so-called elliptic flow 7] vi is a sensitive 
probe of the dense matter produced in a nucleus- nucleus 
collision at ultrarelativistic energies Q- 

While v n , defined by Eq. (JTJ, is a trivial one-particle 
observable which can easily be computed in a model or 
an event generator, the experimental situation is quite 
different. Indeed, the reference direction <j>^ is unknown 
experimentally, and v n can only be measured indirectly, 
from the azimuthal correlations between the detected 
particles. Furthermore, varies randomly from one 
event to the other, which has a remarkable consequence: 



anisotropic flow appears as a truly collective motion, in 
the sense that all outgoing particles in a given event seem 
to be attracted towards some arbitrary direction. 

The standard method for analyzing anisotropic flow is 
to correlate particles with an estimate of &r |9j. How- 
ever, this estimate is itself obtained from the outgoing 
particles , an d one essentially measures a two-particle cor- 
relation Intuitively, two-body correlations are not 
the appropriate tool to probe collective behaviour. In- 
deed, these two-particle methods were shown to be inad- 
equate due to various "nonflow" correlations from quan- 
tum statistics resonance decays, minijet produc- 
tion etc., which are neglected and bias the analysis. 
Recently, new methods were developed, based on higher- 
order (typically, four-particle) correlations, together with 
a cumulant expansion which eliminates low-order nonflow 
correlatio ns 1131 . However, it was argued that experimen- 
tal results could still be biased by nonflow effects 0] 
at this order. In this paper, v n will be analyzed directly 
from the correlation between a large number of particles. 
It will be shown that the results are perfectly stable with 
respect to nonflow correlations, which involve a smaller 
number of particles. 
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Lee- Yang zeroes were already used in analyzing multiplicity dis- 
tributions in high-energy collisions. But it was shown that the 
locations of the zeroes merely reflect general, well-known features 
of these distributions 0, and do not bring any new insight into 
the reaction dynamics. 



FIG. 1: Schematic picture of a nucleus- nucleus collision 
viewed in the plane transverse to the collision axis z. b is 
the impact parameter, $i? its azimuthal angle. cf> is the az- 
imuthal angle of an outgoing particle. 
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Our new method is based on the following global ob- 
servable, which is defined for each event: 



M 

Q = cosn( 



(2) 



where n is the Fourier harmonic under study (n = 1 for 
directed flow vi, n = 2 for elliptic flow), the sum runs 
over all M detected particles, <pj are their azimuthal an- 
gles, and 9 is an arbitrary reference direction. This quan- 
tity is nothing but a projection of the "event flow- vector," 
used in other methods to estimate the orientation of the 
reaction plane , on the transverse direction making an 
angle n9 with respect to the x-axis. In practice, the 
sum in Eq. is often weighted: weights depending on 
the particle mass, transverse momentum and rapidity are 
used in order to reduce statistical errors and increase the 
flow signal. They are omitted here for the sake of sim- 
plicity, but should be included in the actual analysis. 

The central object in the method is the moment gen- 
crating function 
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G(z. 



(3) 



where z is a complex variable, and angular brackets now 
denote an average over a large number of events with the 
same impact parameter. The procedure to obtain v n (as 
will be shown below) is the following: choose a value of 
6; construct Q for each event, evaluate G(ir) for real, 
positive r; plot |G(ir)| as a function of r; determine the 
first minimum tq. The flow is given by v n ~ 2.405/Mro. 

Let us now justify the procedure. We first introduce 
the cumulants Ck, which are defined as [lfij | 



lnG(z) 
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(4) 



The first two terms in this power-series expansion corre- 
spond to the average value of Q , and the square of its 
standard deviation, respectively: 



C i = (Q e ), c 2 = {(Q e f)-{Q e Y 



(5) 



Note that c\ vanishes by symmetry if the detector has 
uniform azimuthal coverage. 

The order of magnitude of the cumulants differs de- 
pending on whether or not there are collective effects 
in the system. Since Q is the sum of M terms of order 
unity, and Ck involves (Q ) k , the naive expectation is that 
Ck should be of order M k or, more generally, scale with 
M like M . As we shall see later, this is precisely the 
case when anisotropic flow is present. When no collective 
effect is present, however, cumulants are much smaller: 
one can view the system as made of independent clusters 
of particles. G(z) then factorizes into the product of the 
contributions of each cluster, which is converted into a 



sum by the logarithm. Hence Ck scales only linearly with 
M. For particles emitted with uncorrelated, randomly 
distributed azimuthal angles, for instance, Eq. JHJ) gives 
c 2 - A//2. 

This shows that the value of Ck for large k is the natu- 
ral observable to characterize collective effects: the larger 
k, the larger the contribution of collective effects, rel- 
ative to other contributions, which scale like M. The 
asymptotic behaviour of Ck for large k therefore provides 
a cleaner separation between collective effects and few- 
body correlations than finite-order cumulants [lflj . This 
can be easily understood physically: the cumulant Ck es- 
sentially isolates the contribution of genuine /c-particle 
correlations, by subtracting out most contributions from 
lower-order correlations. In order to study collective ef- 
fects, which by definition involve a large number of par- 
ticles, k should be as large as possible. 

The asymptotic behaviour of Ck for large k is deter- 
mined by the radius of convergence of the power-series ex- 
pansion, Eq. Q), i.e., by the singularity of lnG(z) which 
lies closest to the origin in the complex plane. Since 
G(z) has no singularity, the only possible singularities of 
lnG(z) are the zeroes of G(z). If z denotes the zero 
closest to the origin, Ck scales typically like Zq for large 
k. Therefore, if Ck scales like M k (collective effects), zq 
scales like 1/M. If there is no collective effect, G(z) is the 
product of contributions of small clusters, and the zeroes 
of G(z) are the zeroes of the individual contributions: zq 
does not depend on M. 

We are now in a position to explain how our approach 
relates to the theory of phase transitions of Yang and 
Lee 0. The starting point is the grand partition func- 
tion: 



+ oo 

00*) = £ Z N e^ kT , 

N=0 



(6) 



where Zn is the canonical partition function for N parti- 
cles at temperature T in a volume V (both T and V are 
fixed). Let /z c denote a reference value of the chemical 
potential //. The probability Pjv to have N particles in 
the system at /i = /i c is 



Pi 



N 



Z Ne ^N/kT 



(7) 



The moment generating function of this probability dis- 
tribution can be simply expressed in terms of the grand 
partition function, Eq. JSJ 



G(z) = PNe zN 

N=0 



Q(H C + kTz) 



(8) 



This function is analogous to our generating function, 
Eq. Js3, with the number of particles N instead of Q , 
and the volume V instead of the multiplicity M. Wc 
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can repeat the previous discussions: if particles are cor- 
related only within small clusters, zq (the zero of G(z) 
closest to the origin) is independent of V. This is the case 
when no phase transition occurs at [i = fi c . Now assume 
that a first-order transition, say, a liquid-gas transition, 
occurs at /i = /i c . Then, the system can be any mix- 
ture of the low-density gas phase and the high-density 
liquid phase. The probability distribution P/v in Eq. (jHJ 
is widely spread between two values N min (gas) and N max 
(liquid) which both scale like the volume V. Then, the 
partition function G{z) depends on the volume V essen- 
tially through the combination zV, and consequently its 
zeroes scale with the volume like 1/V. The general re- 
sult of Lee and Yang is precisely that a phase transition 
occurs at n = // c if the zeroes of G(z) come closer and 
closer to the origin z — as the volume of the system, V, 
increases (note, however, that Ref. is written in terms 
of the variable y — e z instead of z) . 

Let us come back to heavy ion collisions. So far, our 
analysis has been general, and Q could be replaced by 
any extensive variable in Eq. 0. We are now going to 
specify what happens when there is anisotropic flow in 
the system. We can repeat the discussion of Eqs. (|3l5ll . 
but with all average values taken for a fixed orientation of 
the reaction plane <$>r. Such averages will be denoted by 
(• • • |$r). Using the definition of v n , Eq. @J, and sym- 
metry with respect to the reaction plane (which implies 
(smn((/) — 3>h)) = 0), and assuming for simplicity that 
the multiplicity M is the same for all events, one obtains 
from Eq. Q): 



Cl 



(Q e \$ R ) = Mv n cos(n($ R - 0)). 



(9) 



We neglect terms C3 and higher in Eq. J2J. This amounts 
to assuming that the probability distribution of Q is 
gaussian for a fixed this is the central limit theorem, 
which holds if M is large enough, and if there is no other 
collective effect in the system. We further neglect the (In- 
dependence of C2- Finally, averaging over one obtains 
the following theoretical expression of G(z), which we 
denote by G c .\,( z ) since it corresponds to the central limit 
approximation: 



G cX (z) = e c ^^ 2 I (Mv n z), 



(10) 



where Iq is a modified Bessel function. Taking the loga- 
rithm and expanding in powers of z, one checks that the 
cumulant Ck in Eq.Q scales with M like M k , as antici- 
pated. 

The first zeroes of G c ±(z) lie on the imaginary axis at 



ZQ = 17*0 



(11) 



and at — zq, where joi — 2.405 is the first positive root of 
the Bessel function Jq(x). As expected from the general 
discussion above, anisotropic flow v n , being a collective 



effect, is completely determined by z . The situation is 
analogous to a first-order phase transition, in the sense 
that the position of the zero scales like 1/M, and the 
multiplicity M is the analogue of the volume V in Lee- 
Yang's theory. The important difference with statistical 
physics is that the system size is much smaller. As a 
consequence, zeroes never come very close to the origin, 
but the physics involved is essentially the same. 

In a second paper [l7j . Lee and Yang further showed 
that all zeroes lie on the imaginary axis of the variable 
z (or, equivalently, on the unit circle for y — e z ) for a 
general class of models. It is interesting to note that our 
theoretical estimate, Eq. (|lf)|> . has the same property. 
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FIG. 2: Variation of |G(ir)| with r. The crosses are the values 
of \G(ir) | . The solid line displays the expected value \G c .i. (ir) 
defined by Eq. fTTty . 



Figure [21 displays the variation of |G(ir)| as a func- 
tion of r for simulated data. The data set contained 
Aevts = 20000 events. In each event, M = 300 particles 
are emitted independently with an azimuthal distribution 
dN/d(j) oc 1 + 2v 2 cos(2(0 - $r)), where v 2 = 6%, and 
the azimuth of the reaction plane, is randomly cho- 
sen. These numbers are typical values for a mid-central 
Au+Au collision at ^/s^ = 130 GeV, as analyzed by 
the STAR Collaboration ||. The global observable Q 6 
in Eq. Q was constructed for each event with n = 2 
and various values of 9. Figure |21 corresponds to 9 = 0. 
The numerical results are compared with the theoretical 
estimate, Eq. (|10l) . where we have taken c 2 = M/2 (the 
expected value for independent particles). The excellent 
agreement justifies the approximations made in deriving 
Eq. I|10l) . However, a closer look at the numerical results 
(inlay in Fig. EJ) shows that unlike the theoretical esti- 
mate, |G(ir)| does not strictly vanish: due to statistical 
fluctuations, the zeroes of G(z) are slightly off the imag- 
inary axis. This small deviation is physically irrelevant, 
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and we choose to investigate the minima of |G(z)|, rather 
than the zeroes of G(z). We denote by r® the first min- 
imum of |G(ir)|, where the superscript 9 recalls that it 
may depend on the reference angle 9 in Eq. J2J. Identify- 
ing G(z) with the theoretical estimate G c .i.(z), and using 
Eq. (|TT)l . we obtain the following estimate of v n , which 
may also depend on 9: 



3oi 



(12) 
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FIG. 3: The reconstructed value v 2 as a function of 6. 
Squares: simulated data with input v 2 = 6% (same data as 
in Fig. |5J. Circles: simulated data with input V2 = 0. 



This procedure was applied to the simulated data. The 
result is shown in Fig. |3J v 2 coincides with the input 
value v 2 = 6%, up to statistical fluctuations. Performing 
the analysis for several values of 9, and averaging v 2 over 
9, reduces this statistical error. We finally obtain 5.95%: 
v 2 is reconstructed with great accuracy. 

For the sake of illustration, we also applied the same 
procedure to simulated data with no flow. The proce- 
dure yields a spurious flow, due to statistical fluctua- 
tions, which is also shown in Fig. |3J The magnitude of 
this spurious flow can easily be understood. The average 
in Eq. © is evaluated over a finite number of events, 
N cv ts- As a consequence, G(ir) has statistical fluctua- 
tions, whose typical magnitude is l/^/N cvts - For large 
enough 7', they become as large as the expectation value 
given by Eq. i|l(J|) . in which we set v n — and c 2 = M/2. 
This occurs when 

' (13) 



-Mr 2 /A 



As soon as r is larger than this value, fluctuations can 
produce a minimum of |G(ir)|. The corresponding "spu- 
rious flow" given by the analysis, Eq. (|12[1 . satisfies 



< 



V2MlnAT c . 



(14) 



For our simulated data, the right-hand side (rhs) is about 
3.1%, which is depicted as the dashed line in Fig. As 
expected, the values of v 2 lie below this value, but only 
slightly. This is the main limitation of our method: v n 
can be safely reconstructed only if it is larger than the 
rhs of Eq. (QJJI. 

Since no "nonflow" correlation between the particles 
was simulated, standard methods of flow analysis would 
have worked well too. However, the unique feature of the 
present method is its absolute stability with respect to 
such correlations. As an illustration, assume that instead 
of emitting M particles in each event, we emit M clusters, 
each cluster containing q collinear particles. Then, Q e is 
increased by a factor of q. As a consequence, the position 
of the first minimum, 7q , is smaller by a factor of q. Since 
the event multiplicity is now qM, one must replace M 
with qM in the denominator of Eq. 1)12(1 , so that the flow 
estimate v n is strictly the same, as it should. On the 
other hand, estimates of v 2 from 2-particle or 4-particle 
methods |lj| are generally increased by such correlations. 
With the numerical values above, the increase would be 
significant for 2-particle methods (one obtains v 2 = 7.3% 
instead of 6%), but very small for 4-particle cumulants: 
in most cases of interest, these cumulants will give results 
very similar to those obtained with the present method, 
but the latter is the most systematic one to disentangle 
collective motion from other effects. 

The method can be extended to the analysis of dif- 
ferential flow, i.e., the analysis of function of 
transverse momentum and rapidity. This is explained 
in detail in Ref. 0], where we also discuss in detail er- 
rors due to nonflow correlations, statistical fluctuations, 
and show that the method is remarkably insensitive to 
azimuthal asymmetries in the detector acceptance. 

We have shown that Lee- Yang theory of phase tran- 
sitions can be used as a practical means of analyzing 
anisotropic flow experimentally. The method is expected 
to give results similar to cumulant methods, but is sig- 
nificantly simpler to implement, and formally elegant. It 
does not require the knowledge of the reaction plane and 
there is no need to construct correlation functions and 
cumulants. More generally, Lee- Yang zeroes provide a 
natural probe of collective behaviour. It would be inter- 
esting to extend the present approach to other observ- 
ables, in order to look for critical fluctuations which may 
occur in the vicinity of a phase transition. |l9| 



Acknowledgments 

R. S. B. acknowledges the hospitality of the SPhT, 
CEA, Saclay; J.-Y. O. acknowledges the hospitality of 
the Department of Theoretical Physics, TIFR, Mumbai. 
Both acknowledge the financial support from CEFIPRA, 
New Delhi, under its project no. 2104-02. 



5 



[1] C. N. Yang, T. D. Lee, Phys. Rev. 87 (1952) 404. 
[2] C.-N. Chen, C.-K. Hu, F. Y. Wu, Phys. Rev. Lett. 76 
(1996) 169. 

[3] F. Csikor, Z. Fodor, J. Heitger, Phys. Rev. Lett. 82 (1999) 
21. 

[4] Z. Fodor, S. D. Katz, Phys. Lett. B 534 (2002) 87. 

[5] T. C. Brooks, K. L. Kowalski, C. C. Taylor, Phys. Rev. 

D 56 (1997) 5857. 
[6] S. Voloshin, Y. Zhang, Z. Phys. C 70 (1996) 665. 
[7] J.-Y. Ollitrault, Phys. Rev. D 46 (1992) 229. 
[8] STAR Collaboration, K. H. Ackermann, et al., Phys. 

Rev. Lett. 86 (2001) 402. 
[9] P. Danielewicz, G. Odyniec, Phys. Lett. B 157 (1985) 

146; 

A. M. Poskanzer, S. A. Voloshin, Phys. Rev. C 58 (1998) 
1671. 

[10] S. Wang, et al., Phys. Rev. C 44 (1991) 1091. 



[11] P. M. Dinh, N. Borghini, J.-Y. Ollitrault, Phys. Lett. B 
477 (2000) 51; 

N. Borghini, P. M. Dinh, J.-Y. Ollitrault, Phys. Rev. C 

62 (2000) 034902. 

[12] Y. V. Kovchegov, K. L. Tuchin, Nucl. Phys. A 708 (2002) 
413. 

[13] N. Borghini, P. M. Dinh, J.-Y. Ollitrault, Phys. Rev. C 

63 (2001) 054906; Phys. Rev. C 64 (2001) 054901. 

[14] STAR Collaboration, C. Adler, et al., Phys. Rev. C 66 
(2002) 034904. 

[15] Y. V. Kovchegov, K. L. Tuchin, Nucl. Phys. A 717 (2003) 
249. 

[16] N. G. van Kampen, Stochastic Processes in Physics and 
Chemistry (North-Holland, Amsterdam, 1981). 

[17] T. D. Lee, C. N. Yang, Phys. Rev. 87 (1952) 410. 

[18] R. S. Bhalerao, N. Borghini, J.-Y. Ollitrault, Nucl. Phys. 
A 727 (2003) 373. 

[19] S. Jeon, V. Koch, |hep-ph /0304012 



